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Abstract 

In order to test the accuracy of the Heavy Top-mass Expansion (HTE) employed 
in recent two-loop calculations of M w and sin 2 we consider their contributions 
to subtracted quantities of the form {M w ) sub = M W (M H ) - M W (M° H ), where M% 
is a reference point. The results are compared with those obtained by a precise 
numerical evaluation of all the two-loop contributions involving both the Higgs bo- 
son and a fermion loop. For the choice M° = 65 GeV, and over the large range 
65 GeV < M H < 1 TeV, we find very small differences between the precise and HTE 
calculations, amounting to \5M W \ < 0.8 MeV and \Ssl f j\ < 1.2 x 1CT 5 . Although 
corrections involving light fermions are necessary for the consistency and test of ex- 
istent calculations, we also discuss the separate contributions from the top-bottom 
isodoublet. In this case, the differences are larger, although still small, namely 



\SM W \ < 1.9 MeV and \6s 2 eff \ < 4.5 x 10 
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The corrections of 0(g 4 M 2 /M^), evaluated on the basis of Heavy Top-mass 
Expansion (HTE) techniques, are now incorporated in the calculation of the elec- 
troweak observables M w , sin 2 6^, and Tf with f ^ b |]^§. In these papers, the 
HTE is applied in two different ranges of the Higgs boson mass M H , and it is 
found that both expansions match nicely at M H « M w . The results have reduced 
significantly the scheme and scale dependence of the overall corrections, and have 
decreased the estimated Higgs-mass bounds by ~ 30% [Q. 

The dependence of the electroweak observables M w , sin 2 6*^, and on M H 
was also studied recently in Refs. |§, 0. In these calculations, all the two-loop 
contributions involving both a fermion loop and the Higgs boson (H) have been 
taken into account accurately, i.e. without using the HTE, by a combination of 
algebraic and precise numerical methods. For brevity we will refer to this class of 
diagrams as C^ 2 \f, H). We note that graphs of this class in which the H couples to 
external fermions are at most of 0{M 2 /M^ v ) and, therefore, negligible. As not all 
of the two-loop contributions are included, the diagrams in C^(f,H) are actually 
divergent. However, the divergences are M H -independent, so that it is possible to 
evaluate accurately the contributions of (/, H) to subtracted corrections, such as 
(Ar) su b(M H ) = Ar(M H ) — Ar(M°). Here Ar is the correction introduced in Ref. 
and M° is a reference value. C^(f, H) includes all the M H -dependent contributions 
enhanced by factors (M 2 /M^,) 2 (n = 1,2), as well as all remaining M ff -dependent 
effects involving the top quark or a light fermion loop. On the other hand, since two- 
loop purely bosonic self-energy contributions, as well as two-loop boxes and vertex 
parts, are not included, we stress that C^ 2 \f,H) does not contain the full two-loop 
M H - dependence, a limitation that also applies to the calculations of Refs. [JJ-Q]. We 
also point out that, in the approach of Refs. Hf|, the contributions of C^(f,H) 
to subtracted radiative corrections must be calculated at constant values of M w , 
rather than M W (M H ), in order to ensure the cancelation of divergences. The terms 
neglected in this approximation are formally of three-loop order and are expected 
to induce only a small error in the finite parts. 

The aim of this paper is to test the precision of the HTE, as employed in Ref. 
H] U, and the resulting accuracy in the evaluation of C^ 2 \f, H). With this objective 
in mind, we consider the 0(g A ) contributions from Ref. to subtracted corrections, 
relevant to the calculation of M w and sin 2 9^$, and compare the results with those 
obtained for the same quantities from the accurate two-loop calculation of (/, H). 
In order to facilitate the comparison, we employ OSII, one of the on-shell schemes 
of Ref. 0|, since Ref. ||,|6[] employs also the on-shell scheme of renormalization. 
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We also avoid, as much as possible, deviations arising from different treatments of 
higher-order corrections not contained in C^ 2 '(f,H). QCD corrections are excluded 
in both calculations, as they do not play a significant role in the test of the HTE. 
We first note that Ref. || employs the conventional framework || 



2 2 ^ 2 

s 2 c 2 - 



MJ (1 - Ar) 



c 2 = 1 - s 2 = Ml/M 2 z , A 2 = 7ra/V2G 
while Ref. f| uses the alternative expression 

s 2 c 2 = ^(l + Ar N ). (2) 

The subscript N reminds us that, in Eq. (Q), Ar has been introduced in the numer- 
ator. Through two-loop order, Eq. (P leads to 

1 + Ar (l) + ( Ar (l))2 + Ar -(2) ; (3) 



1 - Ar 

where Ar 1 - 1 - 1 is the original one-loop result of Ref. 0, and Af^ stands for a sum 
of explicit two-loop corrections involving the top-bottom isodoublet and the Higgs 
boson (Cf. Eq.(14) of Ref. §). In turn, Ar (1) can be decomposed according to 

Ar (1) = Aa + Ar[l ] + Arj} } + Ar [ x) , (4) 

where Ar^ and ArB denote the one-loop contributions of the top-bottom isodou- 
blet and the light fermions, not contained in Aa, and Ar^\M H ) is the bosonic 
contribution, as defined in Ref. |7j (i.e. including vertex parts and box diagrams). 
The light fermions include the leptons and the first two generations of quarks. At the 
one-loop level the M H dependence resides in Ar^ (M H ). At the two-loop level, the 
M H dependence of Eq. (|) is contained in 2(Ar (1) - Ar[ 1} ) Ar^(M H ) + Af^(M H ) + 
(Ar^\M H )) 2 . The last term, however, does not belong to the C^(f,H) class and, 
for this reason, it is not included in the analysis of Ref. 0. Furthermore, it is not 
affected by the HTE. Therefore, in order to test the HTE by comparing the calcu- 
lations of Refs. H and ||, we disregard (Ar^M^)) 2 . At the two-loop level, the 
contribution of Eq. ([3|) to (Ar7v) s «b(M H ) is then given by 

(Ar N ) sub (M H ) = (Ar^U + (Ar«)L - (Ar«)L + ArfJ,. (5) 
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with the understanding that (Ar b ) 2 sub is not included. The corresponding contri- 
bution in the approach of Ref. || is expressed as 

{Ar N ) sub {M H ) = (1 + 2Aa) (Ar«) su6 + Arfjf + Ar?^, (6) 

where the last two terms stand for the t-b isodoublet and light-fermion contributions 
not contained in Aa. It should be noted that Ar s J b (M H ) in Eq. (|B|) includes all the 
relevant reducible and irreducible two-loop diagrams, while the corresponding light 
fermion contribution in Eq. (||) arises only from the reducible terms in (Ar^) 2 ub . 

In Table 1 we compare the results for (M w ) sub (M H ) obtained by using either 
Eq. (^|) or Eq. (|]), and we list the corresponding shifts 5M W in the case M° =65 GeV 
(the reference point chosen in Ref. |J). We use the input parameters of Ref. [0, 
namely M z = 91.1863 GeV, M t = 175 GeV, (Aa)f = 0.0280. In order to carry out 
the comparison between the two approaches in as close a manner as possible, the 
two-loop contributions in both calculations are evaluated at fixed M w = 80.37 GeV, 
and M W (M H ) is then found by using iteratively Eq. (fj) in conjunction with either 
Eq. (|5[) or Eq. (|6|). (Our conclusions are very insensitive to the precise value of M w 
employed in the evaluation of the two- loop corrections.) The calculation based on 
Eq. (P) employs for M W (M® = 65 GeV) the value obtained from the OSII scheme of 
Ref. ||, subject to the approximation of Eq. (|5|) and the iterative method explained 
above. The quantities (M w ) sub (M H ) (second and third column of Table 1) are ob- 
tained by subtracting M w (65 GeV). It is worth noting that (M w ) sub (M H ) is also 
very insensitive to the precise value of M w (65 GeV). The shift 5M W represents the 
variation of (M w ) sub (M H ) when one employs Eq. @ relative to the value obtained 
from Eq. (|5|). 

From Table 1 we see that, over the large range 65 GeV< M H < 1 TeV, the 
5M W values are very small, |5M W | < 0.8 MeV. There are a number of significant 
differences between the comparison in Table 1 and those carried out in Refs. ||||: 
i) (M w ) sub (M H ), obtained in Refs. ||[| from Eqs.(|,||), is compared with the results 
derived from Eqs.(|2|,[|), rather than Eq. ([[]). In fact, the latter is a resummed expres- 
sion that includes terms of third and higher order involving Aa. The comparison of 
Eqs.(QJ|) and Eqs.(0,|) is much closer, as both expansions are truncated in second 
order and possible deviations arising from different treatments of higher-order cor- 
rections are avoided; ii) As explained before, the contribution (Ar^)^ ufe is excluded 
in Eq. (|5|), in correspondence with Eq. ([|), as it does not belong to C^(f, H) and is 
not relevant to the test of the HTE; iii) The light fermion contribution in Eq. (§) is 
retained, rather than subtracted (the consequence of excluding these contributions 
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in both calculations are discussed later on and in Table 3); iv) As mentioned above, 
in analogy with the treatment of Eq. (||), the two- loop corrections in Eq. (||) are 
evaluated at fixed M w . 

In order to extend these considerations to sin 2 9gft, we recall that, in the on-shell 
renormalization scheme, 



sin 2 6 l e f f (M H ) = k(M H ) s 2 } (7) 



where s 2 = 1 — M^/M 2 , k(M H ) = 1 + Ak is an electroweak form factor, and Ak is 
an important radiative correction. In Eq.(17) of Ref. 0, An is parametrized in the 
form 

8M 2 G„ 

w t \/w «-) -i- - 

B 2' 



Ak(s 2 ) + ^Ap(s 2 ) + AU 2) 



V2 

where the first two terms contain one and two-loop effects, while the third is an ex- 
plicit reducible two-loop contribution. On the other hand, the calculation of Ref. || 
is parametrized in terms of a and s 2 . In the on-shell scheme, physical amplitudes 
are frequently parametrized in terms of and M w (or M z ), as this procedure 
prevents the occurence of large vacuum-polarization contributions involving mass 
singularities 0,0. However, for the purpose of the present comparison, which in- 
volves only subtracted quantities at the two-loop level, it is sufficient to insert in 
Eq. © 8M^Gy V2 = 4vra/[s 2 (l - Ar)], expand (1 - Ar)" 1 « 1 + Ar, and retain 
the additional two-loop contribution (4na/s 2 ) Ar^[Ak^ +c 2 / s 2 Ap^}. In this way, 
the expression based on Ref. is put in a form analogous to that of Ref. , as the 
latter is strictly a two-loop calculation and employs the (a, s 2 ) parametrization. We 
also subtract (47ra/s 2 ) Ar£ [AA^ 1 ) + c 2 /s 2 ApW]{„ since this contribution does not 
belong to the C^ 2 \f,H) class and is not contained in the work of Ref. ||. Numer- 
ically, this term turns out to be very small. As before, in analogy with Refs. 
we evaluate the two-loop contributions at fixed M w . Writing Eq. (0) in the form 
sin 2 O^ff = s 2 + An s 2 , in the approach of Ref. |J the second term is studied con- 
sidering the subtracted quantity (Are s 2 ) su b = (Ak, s 2 )(M h ) — (Ak s 2 )(M°), and 
assuming that (Ak s 2 )(M°) coincides with the value derived from the OSII scheme 
of Ref. f2|. The values of M W (M H ) in each calculation are the ones obtained in the 
analysis leading to Table 1. In particular, s 2 = 1 — M^/Mj in Eq. (0) is evaluated 
in this manner. In Table 2 we compare the subtracted quantity 

(sin 2 6ffiU(M B ) = sin 2 9 l e f f (M H ) - sin 2 9 l e J}(M° H ), (9) 

as evaluated in the two approches for M° = 65 GeV. From Table 2 we see that, over 
the large range 65 GeV< M H <1 TeV, the 5s 2 ^ values are very small, |<5Se//| < 
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1.2 x 10 5 . The results are found to be very insensitive to the precise value of 
sin 2 0$(M°). 

In order to test the accuracy of the approximate treatment of C^ 2 \f,H) in 
Refs. [0-0], which is one of our main objectives, it is important to include in both 
calculations the two-loop effects involving the light fermions. There is also a theoret- 
ical argument that leads to the same conclusion. In a consistent calculation of Ar^ 
at the two-loop level, one should include the reducible contribution 2 Ar^Ar^ in- 
duced when one inserts Eq. (|) into Eq. (§). In fact, without its inclusion, the large 
contribution 2 Act Ar^ , generated by the same substitution, becomes somewhat ar- 
bitrary since, with equal justification, one could separate out Aa + e. For instance, 
— e could be the non-logarithmic part of Aa, or Aa + e could be the MS version of 
Aa, or Aa + e could be the vacuum polarization function evaluated at q 2 ^ Mj. 
With the inclusion of all the fermionic contributions, the ambiguity disappears since 
what is added to Aa must be subtracted from the remaining fermionic component. 
An analogous role is played by Ar^f (M H ) in Eq. @, a contribution that includes 
both reducible and irreducible components. A similar observation applies to the 
calculation of sin 2 O^fl- As explained before, when Ak is parametrized in terms of 
a, reducible contributions proportional to Ar^ are generated, and one must include 
the light fermion contributions in order to obtain an unambiguous answer. How- 
ever, in comparing the two calculations it is also interesting to inquire about the 
specific difference arising from the two treatments of the t — b isodoublet. Indeed, it 
is in these contributions that the enhancement factors (M 2 /M^) n (n = 1, 2) emerge 
at the two-loop level. A simple way of evaluating this difference is to neglect the 
two-loop contributions involving light fermions in both approaches and repeat the 
comparative analysis discussed before. In Table 3 we list the corresponding SM W 
and 5s 2 eff shifts. We find now \5M W \<1.9 MeV and l^s 2 ^ < 4.5 x 1(T 5 , with the 
maximal values attained at large M H . Although these shifts, arising from differences 
in the treatment of the top-bottom isodoublet, are larger than the very small vari- 
ations in the complete calculations, displayed in Tables 1 and 2, they remain small. 
Nonetheless, comparing Tables 1 and 2 with Table 3, we see that, at the O(MeV) 
level in 5M W and O(10 -5 ) in Ss 2 jp the differences in the treatment of the light 
fermions are significant. In fact, their inclusion reduces the magnitude of the shifts 
in the complete calculations. We also stress that, in the on-shell scheme employed in 
this paper, 5M W and 5s 2 ^ are highly correlated. For instance, in Table 3 at M H = 1 
TeV, we have 5M W = 1.9 MeV and this induces a change 5s 2 ^ = —3.7 x 1CT 5 in the 
tree level s 2 = 1 - Ml/M 2 . Thus, the shift 5s 2 eff = -4.5 x 1(T 5 at M H = 1 TeV, 
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shown in Table 3, is mainly due to the effect of 5M W in the tree level correction to 
s^ff, with only a very small change —0.8 x 10~ 5 attributable to differences between 
the precise and HTE evaluation of the radiative correction s 2 Ak. 

As a final check, in order to discriminate the effect of the iteration, we have com- 
pared the calculations of the subtracted radiative corrections (Ar N ) sub and An sub) 
obtained on the basis of Refs. and Ref. 0, when both the one and two-loop 
contributions are evaluated at fixed M w = 80.37 GeV (see also Ref. fllOfl). For 
the differences between the two calculations in the range 65 GeV< M H <1 TeV, 
we find \(ArN) su t)\ < 0.5 x 10~ 4 and |Ak SU 6| < 0.8 x 10 -4 , which induce shifts 
|5M W | < 0.8MeV and |<5s 2 jy| < 1.8 x 10~ 5 , respectively. When the light fermions 
are excluded, we obtain \(Ar N ) sub \ < 1.0 x 10~ 4 and |Ak su;) | < 0.5 x 10~ 4 , which 
induce shifts |5M W | < 1.6MeV and |5Se//l < 1.1 x 10~ 5 , respectively. These effects 
are of the same order of magnitude as shown in Tables 1-3. However, the exact 
details differ, since the fixed M w calculations do not take into account the iterative 
evaluation of M W (M H ) from Eq. (fj) and, for this reason, have less physical meaning. 

From Tables 1 and 2 we see that, over the large range 65 GeV< M H < 1 TeV, 
|5M W | < 0.8MeV, I^Sgjjl < 1.2 x 10~ 5 . The maximal variations are larger, but still 
small, when the light fermion contributions are excluded, namely |5M W | < 1.9MeV, 
\$ s lff\ — 4-5 x 10 -5 . We recall that the current estimate of M H and its 95% C.L. 
upper bounds are in the M H < 300 GeV range. From Tables 2 and 3, we see that in 
that domain the maximum difference in sin 2 9 l ^fj amounts to —1.2 x 10~ 5 (—1.8 x 
10~ 5 ), when the light fermion contributions are included (excluded). Variations of 
this magnitude would induce a change of 2.3% (3.5%), or about 6 GeV (9 GeV), 
in the current 95% C.L. upper bound M H < 262 GeV. On the other hand, the 
latter upper bound [jllj already includes an estimated uncertainty due to higher 
order electroweak effects which is significantly larger than the shifts we have just 
considered. We also note that the choice M° = 65 GeV for the reference point is 
somewhat arbitrary. Although a change in M° does not affect the variation of the 
subtracted quantities between one M H value and another, it does modify |#M W | and 
\Ss 2 eff \ at fixed M H . For instance, if M° = 300 GeV were chosen, 5s 2 e ^ would vary 
from 0.8 x 10~ 5 to 0.5 x 10~ 5 in Table 2, and from 1.1 x 10~ 5 to -2.7 x 10" 5 in 
Table 3, leading to smaller maximal values for |5Sgjj|. 

We conclude our comments with an observation concerning: (Ar^) 2 . Although 
this contribution has not been included in our analysis, it is natural to consider its 
effect in discussions of scheme dependence. For M w = 80.37 GeV, (Ar^) 2 equals 
(0.04,0.11,0.53,1.05,1.57) x 10~ 4 at M H = (65,100,300,600,1000) GeV. At the 
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10 4 level of accuracy, it becomes relevant only at the higher M H values. There 
exists, however, another well-known two-loop contribution involving M H , namely 



the irreducible contributions proportional to [|T^,|T^]. Its effect on Ar is —0.98 x 



10" 4 (M H /TeV) 2 Although it is usually neglected, as in the calculations of 

Ref. ||, our observation is that it would be natural to include it, together with 
(Ar^) 2 , if large values of M H are considered to test the scheme dependence. The 
combined contribution of (Ar^) 2 and the two-loop effects proportional to M 2 equals 
(0.04,0.10,0.44,0.70,0.59) xl0~ 4 . This reduces the magnitude of these effects and the 
ambiguity associated with the possible inclusion or exclusion of their contribution. 
It is also worth noting that, if the M H -dependence of the full two-loop bosonic 
contributions is of the same magnitude as in (Arj, 1 ^) 2 , it would be significantly 
smaller than that arising from the whole C^\f, H). 

In summary, as illustrated in Tables 1-3, by comparing the results of Ref. |2[ with 
those of Refs. |5],|6| in the evaluation of subtracted quantities over the large range 
65 GeV< M H < 1 TeV, we have found only small differences attributable to the use 
of the HTE. When the complete calculations are compared (Tables 1 and 2) they 
are significantly smaller than the errors estimated in Ref. [[|,|J] at fixed M H , while 
they reach about the same magnitude when light fermion contributions are excluded 
(Table 3). We would like to stress, however, that these reassuring conclusions are 
not a substitute for the very difficult, but fundamental task, of achieving a complete 
two-loop evaluation of Ar and other basic radiative corrections of the Standard 
Theory. 
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S 



M H 


(M w ) sub (Eq.5) 


(M w ) sub (Eq.6) 


5M W 


( GeV) 


(MeV) 


(MeV) 


( MeV) 


100 


-23.3 


-22.9 


0.4 


300 


-96.8 


-96.0 


0.8 


600 


-149.8 


-149.7 


0.1 


1000 


-188.1 


-188.8 


-0.7 



Table 1: Comparison of (M w ) m t>(M H ), as obtained from Eq.(|), based on Ref. 
(see text), and Eq.(|5]), based on Ref. [5[]. The latter calculation employs an accurate 
evaluation of the contributing two-loop diagrams, i.e. it does not apply the HTE in the 
two-loop corrections. The shift SM W is the difference between columns 3 and 2. 



M H 


(S 2 eff )sub X 10 3 


(S 2 eff )sub X 10 3 


Ss hf 


( GeV) 


Ref. § 


Ref- IJ 


10" 5 


100 


0.211 


0.207 


-0.4 


300 


0.781 


0.769 


-1.2 


600 


1.152 


1.141 


-1.1 


1000 


1.412 


1.405 


-0.7 



Table 2: Comparison of (sljj-) su b(M H ), as derived from Ref. (see text), with the 
corresponding results from Refs. The latter calculation is based on an accurate 

evaluation of the contributing two-loop diagrams, i.e. does not employ the HTE in the 
two-loop corrections. The shift 8s 2 e ^ is the difference of columns 3 and 2. 
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M H 


5M W 


6s eff 


{M w ) sub 


( s eff)sub 


(GeV) 


(MeV) 


(io- 5 ) 


(MeV) 


(IO" 3 ) 


100 


0.6 


-0.7 


-22.5 


0.202 


300 


1.6 


-1.8 


-94.2 


0.748 


600 


1.9 


-3.9 


-146.5 


1.106 


1000 


1.9 


-4.5 


-184.4 


1.358 



Table 3: The differences SM W and <5Sg// between the calculations based on Refs. 
U and Ref. (see text), when the light-fermion contributions are excluded in both 
analyses. The results reflect the effect of applying the HTE to the two-loop corrections 
involving the top-bottom isodoublet. In columns 4 and 5 we also report the results for 
(M w ) suh and {s 2 e ^) sub from the calculations of Refs. The analogous values in the 

approach of Ref. can be obtained combining columns 2-4 and 3-5. 
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